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We investigate reaction mechanisms occurring in heavy ion collisions at low energy (around 20 
MeV/u). In particular, we focus on the competition between fusion and break-up processes (Deep- 
Inelastic and fragmentation) in semi-peripheral collisions, where the formation of excited systems 
in various conditions of shape and angular momentum is observed. Adopting a Langevin treatment 
for the dynamical evolution of the system configuration, described in terms of shape observables 
such as quadrupole and octupole moments, we derive fusion/fission probabilities, from which one 
can finally evaluate the corresponding fusion and break-up cross sections. The dependence of the 
results on shape, angular momentum and excitation energy is discussed. 



PACS numbers: 25.70.-z, 25.70.Lm, 21.30.Fe, 24.60.Ky 

I. INTRODUCTION 

Nuclear reactions between medium-mass nuclei at low 
energies (around 20 MeV/u) offer the possibility to inves- 
tigate several aspects of dissipative mean-field dynamics 
and to probe nuclear matter under extreme conditions 
with respect to shape, spin and excitation energy. In this 
energy domain, well above the Coulomb barrier but be- 
low the Fermi energies, one essentially observes two types 
of reaction mechanisms: Fusion dominates in the case of 
central and semi-peripheral collisions, while binary re- 
separation processes, associated with deep-enelastic or 
fast fission mechanisms, essentially involve the remaining 
range of (semi-)peripheral reactions flj. However, along 
the transition from fusion to binary processes, composite 
systems with rather elongated shape and large intrin- 
sic angular momentum can be formed, corresponding to 
metastable (or even unstable) conditions, where mean- 
field fluctuations may play a decisive role in determin- 
ing the final outcome. The presence of large event by 
event variances related to the onset of new instabilities 
have been already noted in experiments, from the anoma- 
lous distribution of primary fragment properties in binary 
events 0, The observed variances (in mass, charge, 
excitation energy, angular distribution) appeared much 
larger than the ones predicted by mean-field nucleon ex- 
change models. Similar conclusions were reached in the- 
ory simulations based on stochastic transport models [3] . 

Interaction times are quite long and a large coupling 
among various mean-field modes is expected, leading to a 
co-existence of the different reaction mechanisms in semi- 
central collisions. The study of the competition between 
fusionlike and binarylike processes and, more generally, 
of the fate of the hot nuclear residues created in these 
reactions is a longstanding problem, from which one can 
learn a lot about mean-field dynamics and fundamen- 
tal properties of nuclear forces. This issue has recently 
found a renewed interest, due to the possibility to per- 
form new analyses involving neutron-rich or even exotic 
systems [![. In these conditions the reaction mechanism 



characterizing dissipative collisions is expected to be sen- 
sitive to the density dependence of the isovector part of 
the nuclear interaction, a matter that is largely debated 
nowadays 

In reactions involving medium-heavy nuclei, as a re- 
sult of the complex neck dynamics, one can also observe, 
in sufficiently inelastic collisions, new modes of resepara- 
tion of the colliding system, such as dynamical ternary 
breaking, with massive fragments nearly aligned along a 
common separation axis 0, @ ■ Experimental evidences 
of this mechanism have been recently reported in the 
case of ^^^Au -I- ^^^Au collisions at 15 MeV/u, where 
also aligned quaternary breaking has been observed 
These effects could still be explained in terms of the per- 
sistence of the excitation of shape and rotational modes 
in the projectile- like(PLF) and/or target like(TLF) frag- 
ments that are formed in binarylike events, that would 
lead to further reseparation along a preferential axis, sim- 
ilarly to what happens in fast-fission processes of PLF or 
TLF. It is worth mentioning that, at higher beam energy 
(around 40 MeV/u), where apart from mean- field effects 
two-body correlations are important, ternary breakings 
become the dominant process and new features are ob- 
served, corresponding to the emission of small fragments 
coming directly from the strongly interacting neck region 
[l0| . Actually one may think in terms of a smooth tran- 
sition between the different decay modes of PLF and/or 
TLF, from fast fission, characterized by the splitting into 
fragments with similar size and small relative velocity, to 
neck emission, where small fragments are emitted with 
larger relative velocity with respect to PLF and TLF. 

From the above discussion, it is clear that the un- 
derstanding of the competition between reaction mech- 
anisms in dissipative collisions, as well as of the nature 
of new exotic reseparation modes, requires a thorough 
analysis of the underlying mean-field dynamics and asso- 
ciated shape fiuctuations and rotational effects. In this 
paper, we attempt to improve the dynamical description 
of low energy collisions by coupling a microscopic trans- 
port approach based on mean-field concepts, suitable to 
follow the early stage of the collision up to the formation 
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of composite excited sources, to a more refined treatment 
of the dynamics of shape observables, including the as- 
sociated fluctuations within the Langevin scheme [llj, 
for the following evolution up to the definition of the 
final outcome. In particular, we will discuss the dynam- 
ics of excited sources characterized by given values of 
quadrupole and octupole moments and intrinsic angular 
momentum. This allows one to investigate the compe- 
tition between fusionlike and binarylike reaction mech- 
anisms and to evaluate fusion cross sections, as well as 
the probability and the features of fast-fission processes 
of PLF (or TLF). The paper is organized as it follows: 
In Section 2 we present the hybrid transport treatment 
employed to follow the dynamical evolution of the sys- 
tem. Results concerning the competition between fusion 
and binary processes are discussed in Section 3. Finally 
conclusions and perspectives are drawn in Section 4. 



II. SIMULATION OF THE COLLISIONAL 
DYNAMICS 

A. Dynamical description of nuclear reactions 

The evolution of systems governed by a complex phase 
space can be described by a transport equation, of the 
Boltzmann-Nordheim-Vlasov (BNV) type, with a fluctu- 
ating term, the so-called Boltzmann-Langevin equation 
(BLE) [ll,El: 

I = |[ + {/,i7} = + «[/], (1) 

where /(r,p,i) is the one-body distribution function, or 
Wigner transform of the one-body density, H{r, p, t) the 
mean field Hamiltonian, Icou the two-body collision term 
(that accounts for the residual interaction) incorporating 
the Fermi statistics of the particles, and 6I[f] the fluc- 
tuating part of the collision integral. The nuclear EOS, 
directly linked to the mean- field Hamiltonian iJ, can be 
written as: 

E/A{p, I) = Es/A{p) + Csym{p)l' + 0{I^), (2) 

where / — {N — Z)/A is the asymmetry parameter. We 
adopt a soft isoscalar EOS, Es/A{p), with compressibil- 
ity modulus K = 200 MeV, which is favored e.g. from 
flow studies • For the density (p) dependence of the 
symmetry energy, Csymip), we consider a linear increase 
of the potential part of the symmetry energy with density 
(asystiff): 
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where po is the saturation density, a=13.4 MeV and b=18 
MeV. From the expression of the energy density, Eq.(2), 
the mean-field potential is directly derived. The free 
energy- and angle-dependent nucleon-nucleon cross sec- 
tion is used in the collision integral jl5l |. 



Within such approach, the system is described in terms 
of the one-body distribution function /, but this function 
may experience a stochastic evolution in response to the 
action of the fluctuating term SI[f]. 

However, the numerical resolution of the full BL equa- 
tion is not available yet in 3D. Approximate treatments 
to the BLE have been introduced so far, see Refs.fTH!,!!^, 
such as the Stochastic Mean Field (SMF) model, that 
consists in the implementation of stochastic density fluc- 
tuations only in coordinate space and can be solved nu- 
merically using the test particle method [15j. The lat- 
ter approach has shown to be particularly appropriate 
for the description of the evolution of the dilute unsta- 
ble sources that develop in dissipative collisions at Fermi 
energies (30-100 MeV/u) [H. However, here we are 
essentially interested in semi-central reactions at lower 
energies where, most likely, the formation of elongated 
(rather than dilute) systems is observed, and phenomena 
associated with surface (rather than volume) metastabity 
and/or instability may take place. To improve the treat- 
ment of fluctuations suitable to describe the latter sce- 
nario, we will adopt a hybrid description of the dynamics: 
We follow the miscroscopic SMF evolution until the time 
instant when local thermal equilibrium is established and 
one observes the formation of quasi-stationary elongated 
systems, with density close to the normal value. Then, to 
deal with the following evolution of the system, we move 
to a more macroscopic model description, where the sys- 
tem is characterized in terms of global observables, for 
which the full treatment of fluctuations in phase space is 
numerically affordable, as explained below. 



B. Dynamical evolution of shape observables 

This Section is devoted to the description of the dy- 
namical evolution of excited systems whose leading de- 
grees of freedom are shape observables, while the density 
keeps always close to the normal value, po = , being 

To the nuclear radius constant (rg = 1.2 fm). The con- 
figuration of the system under study, having given charge 
Z and mass A, is described by three global observables 
(and associated velocities): the quadrupole moment /32, 
the octupole moment and the rotation angle uj. For 
situations far from the spherical shape the thermal agi- 
tation can induce fluctuations that may eventually lead 
to break-up channels. Hence the correct treatment of 
shape fluctuations is crucial for the characterization of 
the reaction mechanism. To this purpose, we consider 
the stochastic extension of the Ray leigh- Lagrange equa- 
tions of motion fT§| (the Langevin equation): 
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where qi(i — 1,2,3) 

Ekin{(liAi) + Erot{<liAi 



(w,^2,/?3)- L{liAi) = 
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grangian of the system and 
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is the Rayleigh dissipation function. Ekin , Erot and Epd 
indicate the kinetic, rotational and potential energy of 
the system, respectively, and the quantity Rij is the dis- 
sipation tensor. The difference with respect to the stan- 
dard Rayleigh-Lagrange equations is the fluctuation term 
Ffiuc, that can be interpreted as a rapidly fluctuating 
stochastic force, in the same spirit of the Brownian mo- 
tion, similar to the fluctuating term of the BLE, Eq.(l). 
We solve numerically the set ^ of coupled equations. 

For given values of the quadrupole and octupole mo- 
ments, the shape of the system is parametrized, in terms 
of the polar angle 9, as it follows: 



R{0) = Ra{f32,P3){l + f3i{f32,f33)Yio{ey 



(6) 



where the functions l^io(^) are spherical harmonics. The 
parameters /3i and Rq are introduced to conserve the 
position of the center of mass and the total volume V of 
the system and can be determined from the equations: 
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T 

dV =— R^{9)sin 
-J Jo 



i?'*(6') sin 61 cos 6'd6' = 0, (7) 



(8) 



where z denotes the coordinate along the system maxi- 
mum elongation axis (or symmetry axis) . In the following 
we discuss in detail the derivation of the different terms 
of the Lagrangian L. 



To calculate the mass tensor Mij , we adopt the prescrip- 
tions of Ref. 191: 



(12) 



with 



M;^. =2^mpo / y^huR'+^{9)Pi{coB9){'^S + 
Jo oPj 

+i?o (^|^Yio + l",o)}sin0d0.(13) 

Here Pi are Legendre polynomials and S{9) = In 
our calculations we have L = 5. The coefficients b2i and 
63/ are obtained solving the system of equations: 

L 

Xl^fe'^™' = k^l...L, m = 2,3 (14) 
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with 



^ki= R'-\9) l^lPi{cos9) 



1 dR{9)dPi{cost 
'R(9)~d9 d9~ 



(15) 



R''-\9)l^kPk{cos( 



1 dR{9) dPk{cost 
'R(9)~d9 dr~ 



1. Rotational energy 
The rotational energy is simply equal to: 

Erot = llW2,P3)LO^ (9) 

where 

/(/32,/33) = r R''{9){l+cosH}sm9d9 (10) 

5 Jo 

is the moment of inertia for the whole system, being m 
the nucleon mass. 

2. Kinetic energy 

The kinetic energy can be expressed as it follows: 
1 ^ 

Emu = 2 E M,j{l32,P3)qtqj (H) 



3. Potential energy: Nuclear term 

Concerning the nuclear part of the potential energy. En, 
we discuss essentially the surface contribution, since our 
system keeps a volume constant in time. We adopt a 
double volume integral of the Yukawa-plus-exponential 
folding function (20| : 

aJl - k,P) f f fa \ e-"!'' „ „ , 

En^~ 'I 2 2 3 W / --2 dVdV, 

(17) 

where is the surface-energy constant, kg is the surface- 
asymmetry constant and a is the range of the Yukawa- 
plus-exponential potential, a denotes the modulus of the 
relative distance a = \r — r'\. Parameters have been fitted 
to the ground-state energies and fission barrier heights 
(21I [2^ . In order to reduce the numerical efforts, the 
integral of Eq. p7| can be transformed into a double sur- 
face integral, by using the twofold Gauss divergence theo- 
rem. For axially symmetric shapes, one of the azimuthal 
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integrations can be performed trivially [20|, |23| and the 
resulting threefold integral is: 



as{l - kgP 



2 cr 

+ 2- + 2 
a 



-a I a 



5. Dissipation function 

The one-body dissipation mechanism is evaluated as it 
follows (see Ref.[T§| for details): 



dE 
dF 



mpov (p fi dS 



(21) 



(^ l^^ where the integration is performed over the whole surface 
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of the system, v — jVp is the average nucleon velocity 



and 



tensor Ri 



where the distance a can be expressed as: 

<T = [i?2(6») + r2(9') - 2R{e)R{e') 

■ {cos 9 cos 9' + sin 9 sin 9' cos f/'}] 

and 

P{9,9',(P) = R{9)sm0{R^{0)-R{9)R{0')[cos9cos0'- 



sin0sin9' cosS-R{9')^^^\sm9cos9'-cos0siTi9' cosS} ^ ^ 
d9 



\d]Z\ 
I dt I 



r, ■R^r-R{9) 



(22) 



Hence we get the following expressions for the dissipation 



Ri^j = 'Z'' ■"nmpov 
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dRo 

dp, 



1 f dR{e) \ 

W)\ — 



d9 J 

S\ R^{9)sin0d0 (23) 



4.. Potential energy: Coulomb term 
The Coulomb part of the potential energy is taken as 



6. The Langevin term 

The stochastic force Ffiuc{t) will determine fluctua- 
tions in momentum space, according to the value of the 
diffusion coefficient D. We assume that 
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C 
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where jt;^°-^p jg the Coulomb energy corresponding to a 
sharp charge density distribution and AS* is a correc- 
tion due to the diffuseness. 

The sharp-surface part of the Coulomb energy is equal 
to: 



E. 



sharp Pp^ 

^ 6~ 



P{9,0',^)P{0',0,-^) 



d0d0'f 



(19) 

where pp is the charge (proton) density, pp — ^po- The 
correction to the Coulomb energy due to the diffuseness 
can be expressed as: 
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where ac is the range parameter of the Yukawa function 
generating the diffuse charge distribution p3l - l25| . 



(24) 



The action of the stochastic force Fji^c may be simulated 
numerically by repeatedly producing a random kick 6P 
in the collective velocity associated with the quadrupole 
and octupole moments. The value of bP is chosen ran- 
domly from a Gaussian distribution with mean value and 
variance given by: 



,5P = 



((5P)2 = D5t 



(25) 



(26) 



where 8t is the small time step between two kicks. The 
diffusion coefficient D can be found using the Einstein 
relation: 



D = 2T7, 



(27) 



where 7 is the dissipation coefficient and T is the tem- 
perature of the system [2^ . Hence the fluctuations that 
we are considering are induced essentially by the ther- 
mal agitation. We notice that our dissipation tensor Rij , 
introduced above, has also nondiagonal terms. Hence, 
to correctly extract the dissipation coefficients, we diag- 
onalaze the dissipation tensor Rij 7^^ . The tensor 7^^ 
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will have only diagonal elements: 72 and 73. Now we 
can find D2 and in the new coordinate system and 
evaluate 5P2 and 5P^, the random kicks for the new co- 
ordinates. Finally it is possible to go back to the general 
coordinates /32 and /^a, by the inverse transformation, and 
obtain SPp^ and SPp^. 



III. RESULTS 

We will exploit the Langevin treatment outlined above 
to investigate the competition between (incomplete) fu- 
sion and binary break-up mechanisms in low energy reac- 
tions. We consider the system ^®Ar -I- ^^Zr at two beam 
energies, 9 and 16 MeV/u, in the following range of im- 
pact parameters: b = 5-7 fm and b = 4-6 fm at 9 and 
16 MeV/u, respectively. Within this selection, accord- 
ing to the SMF dynamical evolution, one observes the 
formation of rather elongated configurations for which 
fluctuations are expected to be crucial in determining 
the following evolution. For lower impact parameters, 
the conditions of the reactions are such that one always 
obtains incomplete fusion, while for larger impact param- 
eters bynary break-up is observed. Contour plots of the 
density in the reaction plane, as obtained in the SMF 
calculations, are displayed in Figs. 1-2, for the two reac- 
tions. 

The description of the system in terms of the global 
observables /32, Pi and uj begins at the moment when, 
according to the full SMF evolution, the composite sys- 
tem reaches a quasi-stationary shape, having dissipated 
almost completely the radial part of the kinetic energy 
deposited into the system, while the angular part is con- 
verted into intrinsic spin. This time instant is estimated 
to be around t freeze-out ~ 200 fm/c. During the earlier 
dynamical evolution, pre-equilibrium nucleon emission 
takes place. As a consequence, mass and charge of the 
system are smaller than the total mass and charge num- 
bers, respectively. We get A w 122, Z w 53. As one can 
see from Figs. 1-2, the system configuration can be suit- 
ably parametrized in terms of quadrupole and octupole 
moments. From this point of view, the Langevin treat- 
ment introduced above appears appropriate to describe 
the following evolution, though the dynamical descrip- 
tion is devolved to few leading degrees of freedom. The 
initial conditions of the Langevin equation have been de- 
termined running 10 SMF trajectories. The correspond- 
ing parameters are listed in Tables LII, for a couple of 
events, for each considered case. 

Then, within the Langevin treatment, 200 stochastic 
events were considered for each SMF trajectory. Fluctu- 
ations are injected each 3 fm/c. 

According to the values listed in Tables LII, we test 
essentially the behavior of composite systems with a va- 
riety of conditions of angular momentum, ranging from 
50 H to 100 h and quadrupole moment (32, from 0.2 to 
1. The excitation energy is about 250 MeV, correspond- 
ing to temperatures of the order of 4 MeV. Apart from 



b (fm) 


E* (MeV) 


L (h) 


P2 


/33 


dl32 
dt 


dt 


P 


7 


225 


100 


1.14 


-0.73 


0.099 


0.024 


0.990 


7 


242 


95 


1.00 


-0.76 


0.143 


-0.129 


0.990 


6 


240 


77 


0.83 


0.47 


0.062 


-0.010 


0.645 


6 


224 


84 


1.01 


-0.52 


0.113 


-0.063 


0.880 


5 


216 


64 


0.58 


-0.32 


0.125 


0.938 


0.375 


5 


227 


58 


0.56 


0.36 


-0.004 


0.005 


0.145 



TABLE L Characteristics of the composite system, as ob- 
tained in the reaction ^"Ar + ^''Zi at 9 MeV/u at the time 
t freeze-out- Excitatlon energy, intrinsic angular momentum, 
quadrupole moment, octupole moment and associated collec- 
tive velocities. The time unit adopted to define the collective 
velocities is 10~'^^s = 30 fm/c. Two events are displayed for 
each impact parameter. The fission probability (see text) is 
reported in the last column. 



b (fm) 


E* (MeV) 


L (ft) 


P2 




d/32 
dt 


d/33 
dt 


P 


6 


279 


90 


0.88 


0.34 


0.016 


0.059 


1.000 


6 


277 


97 


0.88 


0.44 


-0.015 


-0.031 


1.000 


5 


241 


73 


0.37 


0.15 


-0.063 


-0.047 


0.320 


5 


252 


77 


0.63 


0.40 


0.136 


-0.020 


0.580 


4 


258 


63 


0.31 


0.06 


0.052 


0.018 


0.110 


4 


247 


52 


0.22 


0.05 


-0.007 


0.002 


0.035 



TABLE II: Same as in Table I, but for the reaction at 16 
MeV/u. 



the situation observed in the case of b = 7 fm, E/A = 
9 MeV/u, the octupole moment, /Ss, always takes rather 
small values, of both signs, indicating that the memory 
of the entrance channel mass asymmetry is lost. Also the 
quadrupole and octupole collective velocities are rather 
small and may take values of both signs, suggesting that 
collective motions, apart from the rotation associated 
with the intrinsic spin, are damped. These conditions 
correspond closely to quasi-stationary, metastable situa- 
tions, i.e. the system is stable against small shape fluc- 
tuations. From one side, it may evolve radiating its ex- 
citation energy and spin and relaxing slowly towards the 
spherical configuration. On the other hand, if the ampli- 
tude of the kicks of the associated collective velocities is 
large enough, the system may overcome the fission bar- 
rier and reach configurations corresponding to surface in- 
stabilities, from which it rapidly separates in two pieces. 
However, one should also consider that the latter possibil- 
ity is in competition with nucleon emission, that reduces 
the excitation energy (and the associated amplitude of 
thermal fiuctuations), while the shape of the system is 
evolving. The nucleon emission rate can be evaluated 
according to the standard Weisskopf formalism [23] . For 
the situations under study, the excitation energy reduces, 
due to nucleon emission, approximately by 2.5 MeV each 
30 fm/c. We follow the trajectory of the system until the 
available excitation energy is fully dissipated. 

Hence, thanks to the introduction of fluctuations in 
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FIG. 1: Contour plots of the density projected on the reaction plane calculated with SMF for the reaction Ar + Zr at 9 
MeV/u, at several times (fm/c). The size of each box is 40 fm. 
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FIG. 2: The same as in Fig.[T]but at 16 MeV/u. 



the dynamical evolution, for a given impact parameter 
one observes a bifurcation of trajectories, leading either 
to compact shapes (fusion) or to elongated shapes, with 
large values of quadrupole and/or octupole moments, 
that eventually cause the break-up of the system. Actu- 
ally the two possible outcomes are associated with a kind 
of bimodal behavior of the shape observables, related to 
configurations corresponding to local minima of the to- 
tal (surface -I- Coulomb) energy. In may be interesting to 
notice that bimodality has been recently observed also in 
the context of liquid-gas phase transitions, where volume 
instabilities are concerned and dilute systems may either 
recompact to normal density or split into a huge number 
of small fragments [28!|. 



A. Fission rates 



In the following, we will first discuss some illustrative 
results obtained in the case of the reaction at 16 MeV/u, 
b — 5-6 fm. In Fig. 3 we present one example of trajecto- 
ries corresponding to the two possibile exit channels (fu- 
sion or fission), in the (/32, /^s) plane. Due to the random 
kicks, starting from the same initial conditions, rather 
different paths are explored. It should be noticed that, 
also in the case of trajectories leading to fusion, the final 
configuration is not exactly spherical, but is associated 
with small (not vanishing) values of the quadrupole mo- 
ment. This corresponds to the stationary configuration 
compatible with the amount of intrinsic angular momen- 
tum present in the system. On the other hand, break-up 
configurations are characterized by rather large values of 
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FIG. 3: One example of trajectories leading either to fusion or 
to break-up, in the {132, Pi) plane, as obtained in the reaction 
at 16 MeV/u, b = 5 fm. 
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FIG. 4; Correlations between the values of quadrupole and 
octupole moments, as obtained for the break-up configura- 
tions in the case of the reaction at 16 MeV/u, b = 5-6 fm. 



/32 and/or fi^. Actually one sees an interesting correlation 
between the two parameters, that is represented in Fig. 4. 
In fact, both a large quadrupole or octupole moments 
are linked to break-up configurations, that correspond to 
tangent spheroids. Fluctuations of the octupole moment 
are rather large, though the majority of the events is lo- 
cated near = 0, corresponding to symmetric fission. 

In Figs. 5-6 (left) the fission rate, dN/ dt, as obtained for 
& = 6 - 7 /to at 9 MeV/u and 6 = 5 - 6/to at 16 M eV/u, 
is displayed as a function of time for a set of 200 events 
in each of the cases considered. For the most peripheral 
impact parameters, after an initially increasing trend, re- 



lated to the time interval needed to build and propagate 
fluctuations, we observe an almost exponential decrease, 
as expected in the case of constant break-up probability 
Ibreak- In this case one can write: dN/dt = Ntjbreak, 
with Nt = Aroe-T'"-"'=* and A^o = 200 (the total num- 
ber of events considered) . This corresponds to situations 
where the break up probability {"Streak ~ 0.002 c/ fm) 
is not much affected by the competing nucleon emission. 
AU events practically lead to fission over a time interval 
that is shorter than the one needed to exhaust the avail- 
able excitation energy by nucleon emission. In fact, the 
maximum of the emission rate is observed at about 300 
fm/c and the system needs, on average, roughly 500 fm/c 
to reach the break-up configuration (this is actually the 
half life time Ttreak = l/76reafc)- On the other hand, for 
smaller impact parameters (corresponding to lower defor- 
mation of the system and lower angular momentum) , the 
break-up probability jbreak is quenched approximately 
by a factor 4 (see the left panel of Figs. 5-6) and de- 
creases in the course of time because of nucleon emission, 
that reduces the excitation energy and the correspond- 
ing amount of thermal fluctuations. In most cases, the 
excitation energy deposited into the system is dissipated 
before the break up conflguration may be reached. It is 
interesting to notice that, even in the most favourable 
case, the typical times of the process are rather large 
(500 fm/c), as compared for instance to the time scales 
associated with the development of volume instabilities 
in multifragmentation processes at higher energies (about 
150 fm/c). This can be explained in terms of the larger 
amount of excitation energy deposited into the system in 
the latter case (that induces fluctuations of higher ampli- 
tude and collective radial expansion) and of the smaller 
growth times associated with volume instabilities [l3| . 

The corresponding fraction of events that undergo 
break-up, Phreak, is reported in Table I-II, at the two 
energies and for all impact parameters considered. From 
the estimated break-up probabilities it is possible to con- 
struct the fusion cross-section, (7/(6) = [l—Pbreak) ^nbdb, 
that is displayed in Fig. 7, for the two energies. We also 
show, for comparison, the results obtained within the 
SMF approach only, where due to the approximate treat- 
ment of fluctuations, one gets distributions close to a 
sharp cut-off (approximated by a sharp cut-off in the flg- 
ure). It is interesting to notice that, especially in the case 
of the reaction at 9 MeV/u, the fusion cross section is 
reduced significantly by the introduction of fluctuations 
that, in turn, help the system to overcome the fission 
barrier and to break-up. 



B. Features of fission fragments 

The time ttreak, needed to reach the break-up config- 
uration, is connected to other interesting features of the 
reaction dynamics, depending on the various entrance 
channel conditions. In fact, due to the intrinsic spin, 
the system rotates while its shape evolves according to 
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FIG. 5: (Color online) Left panel: Distribution of the time 
tbreak (§66 text), as obtained for the reaction at 9 MeV/u and 
impact parameters b = 7 fm (full histogram) and b = 6 fm 
(shaded histogram). Right panel: Angular distribution of the 
break-up direction. Full line and circles refer to b = 7 fm; 
dashed line and open squares are for b = 6 fm. 



to the rotation (on the plane perpendicular to the direc- 
tion of the intrinsic spin of the system) until the break- 
up configuration is reached is shown in Figs. 5-6 (right 
panel), for the two energies and two impact parameters. 
Obviously, the shape of this distribution depends on the 
fission probability, but also on the system angular veloc- 
ity (that in turn depends on the intrinsic spin). In fact, in 
absence of rotation (vanishing spin) the fragments would 
always be emitted along a fixed axis. In the case of the 
most peripheral events, a clear peak is observed in the 
distribution. On the other hand, for more central impact 
parameters, the half life time is much larger and one es- 
sentially gets a flat distribution for Otreak, similarly to 
what is expected in the case of standard statistical fis- 
sion. 



C. Fast-fission of PLF and TLF 
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Several shape, angular momentum and excitation en- 
ergy conditions can be observed also, in the case of col- 
lisions between heavy systems, after the separation into 
PLF and TLF, for one (or both) of these products. Thus 
it is interesting to investigate fast fission processes of 
these objects, leading to ternary (or quaternary) breaking 
of the whole system. For instance, we display in Fig. 8, 
density contour plots as obtained in SMF simulations of 
semi-peripheral collisions of An + Au at 15 MeV/u, for 
which aligned ternary and quaternary breaking has been 
recently observed experimentally Q. One can see that 
similar shape configurations, as the ones observed in the 
reactions investigated above, may appear for PLF/TLF 
fragments. However, these fragments have lower angu- 



FIG. 6: (Color online) Same as Fig.5, for the reaction at 16 
MeV/u and b = 5-6 fm. 



Eq.(4). As a consequence, the direction along which 
the system separates into pieces is strictly connected to 
threak- Hcncc the shape of the angular distribution of 
fission fragments can be used as a clock of the collision, 
from which one can extract information on the break- 
up probability and the underlying reaction mechanism. 
This is an appealing issue that can be investigated also 
experimentally by looking at the angular distribution of 
the emerging reaction products and at the possible ex- 
istence of alignment effects [1, [l3| • In the case of a fast 
break-up (fast fission) the angular distribution should ex- 
hibit a peak: due to the elongated shape of the system, 
the emission is not isotropic. Along the separation pro- 
cess, fragments acquire velocities essentially due to the 
Coulomb repulsion, according to the Viola systematics, 
like in standard fission, but with a preferential emission 
axis. The distribution of the angle, Ohreak, corresponding 
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FIG. 7: (Color online) Fusion cross section, as a function 
of the impact parameter b, as obtained in the reactions at 
9 MeV/u (black histogram) and 16 MeV/u (grey histogram) 
with the Langevin treatment, Eq.(4). The lines correspond to 
SMF simulations at 9 MeV/u (full) and 16 MeV/u (dashed). 
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FIG. 8: Contour plots of the density projected on the reaction plane calculated with SMF for the reaction ^^^Au + ^''^Au at 
15 MeV/u, b = 6 fm, at several times (fm/c). The size of each box is 80 fm. 



lar momentum (about 20 — AOh) and excitation energy 
(of the order of 100 MeV). The corresponding break-up 
probabihty is of the order of 10% and emission times 
are longer (« 2000 fm/c). The fast-fission mechanism 
could explain qualitatively some of the features observed 
experimentally, such as alignment effects and fragment 
relative velocities and charge distributions. However, a 
thourough analysis of the kinematical properties of the 
reaction products [9], as well as the estimated rather 
short break-up times, suggests the persistence of non- 
equilibrium effects in momentum space, i.e. the presence 
of collective velocities in (32 and/or fS^, in addition to the 
tangential velocity generated by the intrinsic angular mo- 
mentum. Collective velocities, probably underestimated 
in the SMF calculations, would speed up the fragmen- 
tation process since the system is pushed towards more 
exotic shapes, from which it is easier to overcome the 
fission barrier. 



IV. CONCLUSIONS 

In this article we have investigated the role of shape 
fiuctuations in the dynamical evolution of excited sys- 
tems that can be formed in semi-peripheral reactions at 
low energy (around 20 McV/u). Quasi-stationary com- 
posite systems, with quadrupole and/or octupole defor- 
mation, are observed, for which shape fluctuations are 
essential to overcome the fission barrier and eventually 
break-up. This analysis is performed within a hybrid 
treatment that couples the study of the early stage of the 
dynamics, devolved to a microscopic stochastic transport 
approach, up to the formation of primary excited sources, 
to a full Langevin description of the leading degrees of 



freedom of these objects: quadrupole, octupole moments 
and angular velocity. For temperature, shape and an- 
gular momentum conditions obtained in semi-peripheral 
reactions, typical time scales of the break-up process are 
of the order of 500 fm/c. The fission fragments are emit- 
ted along a preferential direction, that corresponds to 
the maximum elongation axis. Due to angular momen- 
tum effects, this direction may rotate while the shape of 
the system is evolving towards break-up configurations. 
Hence a careful analysis of the angular distribution of the 
reaction products may give relavant information about 
fission probabilities and the involved time scales, that in 
turn are closely linked to the mean-field dynamics and 
the properties of the nuclear interaction (range, surface 
energy, two-body correlations) . From this study it is clear 
that a good treatment of mean-field fluctuations is a cru- 
cial point in the characterization of dissipative reactions. 
The model employed here provides a suitable description 
of surface modes, parametrized in terms of quadrupole 
and octupole oscillations, but it could miss some non- 
equilibrium effects that can help the system to break- 
up. In fact, collective velocitities related to shape ob- 
servables are likely underestimated in the SMF approach 
[29| and the role of multipolarities higher than octupole 
is neglected in the Langevin treatment. A fully micro- 
scopic description of the whole process would be highly 
desirable, though it is far from being trivial. Some at- 
tempts are represented by improuved quantum molecular 
dynamics calculations (ImQMD) (sO]. Stochastic exten- 
sions of Time-Dependent-Hartree-Fock (TDHF) calcula- 
tions should also provide a valuable tool to characterize 
reaction mechanisms in low energy collisions [sij . Work 
is in progress in this direction. 
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